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On the continuum limit of the Discrete Regge model in 4d 

Elmar Bittner''''^'Q, Wolfhard Janke'^'S and Harald Markum'^ 
'^Atominstitut der Osterreichischen Universitaten, A-1040 Vienna, Austria 
''Institut fiir Theoretische Physik, Universitat Leipzig, D-04109 Leipzig, Germany 

The Regge Calculus approximates a continuous manifold by a simplicial lattice, keeping the connectivities of the 
underlying lattice fixed and taking the edge lengths as degrees of freedom. The Discrete Regge model employed 
in this work limits the choice of the link lengths to a finite number. This makes the computational evaluation of 
the path integral much faster. A main concern in lattice field theories is the existence of a continuum limit which 
requires the existence of a continuous phase transition. The recently conjectured second-order transition of the 
four-dimensional Regge skeleton at negative gravity coupling could be such a candidate. We examine this regime 
with Monte Carlo simulations and critically discuss its behavior. 
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1. REGGE QUANTUM GRAVITY 

One promising method to quantize the theory 
of gravitation employs the Euclidean path inte- 
gral 



(1) 



where the partition function describes a fluctu- 
ating space-time manifold. In the Regge ap- 
proach the (quadratic) link lengths qi represent 
the dynamical degrees of freedom deforming 
continuously a simplicial lattice with fixed inci- 
dence matrix, whereas in the somehow comple- 
mentary approach of dynamical triangulation the 
incidence matrix is fluctuating with constant link 
lengths. Any of these two approaches is plagued 
with various problems; for the Regge approach 
see[|§. 

In "conventional" Regge theory the Regge- 
Einstcin action including a cosmological term , 



(2) 



is used. The first sum runs over all products of 
triangle area At times corresponding deficit angle 
weighted by the bare gravitational coupling /?. 
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The second sum extends over the volumes of 
the 4-simplices of the lattice and allows together 
with the cosmological constant A to set an overall 
scale in the action. 

The Discrete Regge model was invented in an 
attempt to reformulate (ih as the partition func- 



tion of a spin system 



It is defined by re- 



stricting the squared link lengths to take on only 
two values 



(3) 



where hi = 1,2,3, and 4 for edges, face diago- 
nals, body diagonals and the hyperbody diagonal 
of a hypercube, respectively, is chosen to allow 
for fluctuations around flat space. The Euclidean 
triangle inequalities are fulfilled automatically as 
long as e is smaller than a maximum value Cmax- 
The measure D[q] in the quantum gravity path 
integral is taken to unity for all possible link con- 
figurations. 

Numerical simulations of the Z2 system become 
extremely efficient by implementing look-up ta- 
bles and a heat-bath algorithm. In this work typ- 
ically 200 000 — 500 000 iterations have been gen- 
erated. Calculations have been performed with 
the parameter e — 0.0875 and the cosmological 
constant A = because (||) already fixes the av- 
erage lattice volume. 
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Figure 1. FSS of the specific- heat maxima C„ 
as a function of the lattice size L. 
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Figure 2. FSS of the susceptibihty maxima Xmax, 
indicative for a cross-over. 



2. RESULTS 

The Discrete Regge model - like full Regge the- 
ory - exhibits two phase transitions Q. One is 
located at a negative value and the other one at a 
positive value of the bare gravitational coupling. 
Although earlier work concentrated on the latter 
transition, there is no reason for favoring a posi- 
tive value of (3 over a negative one. Since a Wick 
rotation from the Lorentzian to the Euclidean sec- 
tor of quantum gravity is not feasible in general, 
the sign in the exponential of the path integral is 
not fixed a priori. 

The transition for positive gravitational cou- 
pling has been observed before the one at negative 
/3 1^. There, histograms, e.g. of AtSt @], clearly 
show a two-peak structure. The two phases coex- 
ist and tunneling from one phase to the other and 
back takes place. The system also exhibits a hys- 
teresis; the transition occurs at a larger value of 
(3 if the simulation is started from a configuration 
in the well-defined phase than it does if the cal- 
culation is started from a "frozen" configuration. 
Given all these pieces of evidence, we conjecture 
that the transition at /3 > is of first order. 

To determine the order of the transition at 
negative f3 we employed in Ref. histogram 
techniques and used the Binder-Challa-Landau 
(BCL) cumulant criterion jo). The BCL cumu- 

lant is defined as Bl := 1 — ^^^2^2 ; with E being 
the action of the system under consideration. It 
was evaluated for AtSt on different lattice sizes 



with L = 3 to 10 vertices per direction, simulat- 
ing the system at several values of the bare cou- 
pling jS with high statistics. For the BCL cumu- 
lant a trend towards 2/3 was observed and all his- 
tograms showed a clear one-peak structure, cf. 0. 

In the present Monte Carlo simulations we 
recorded for every run the time series of the 
energy density e = E/Nq and the magnetiza- 
tion density m = ai/Ng, with the lattice size 
A'o — L^- To obtain results for the various ob- 
servables O at values of the bare gravitational 
coupling P in an interval around the simulation 
point /3o, we applied the reweighting method [T^ . 

With the help of the time series we can compute 
the specific heat, 

C{(3)^rj'No{{e')~{e)') , (4) 

and the (finite lattice) susceptibility, 

X{f3) = Noiim') ~ {\m\f) . (5) 

Figures |l| and || show the finite-size scaling (FSS) 
of the maxima of the specific heat Cmax and the 
susceptibility Xmaxj respectively. While the be- 
havior of Cinax could Still be explained as critical 
scaling with a negative exponent a, the flatten- 
ing of Xmax is rather unusual for a second-order 
transition and should be taken as an indication 
for a cross-over regime. 

Another feature of the system can be seen 
in Figs. H and ^ depicting the critical gravita- 
tional coupling, as determined from the speciflc- 
heat maxima Cmax and the susceptibility maxima 
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Figure 3. Volume dependence of the critical grav- 
itational coupling /3(Cmax), as determined from 
the specific-heat maxima Cmax in Fig. |^. 




Figure 4. Volume dependence of the critical grav- 
itational coupling /3(xmax), as determined from 
the susceptibility maxima Xmax in Fig. ^. 



Xmax- In the case of a second-order transition the 
extrapolations of all pseudo-transition points lead 
to one infinite- volume critical value. This seems 
to be violated in the four-dimensional Discrete 
Regge model and thus favors the interpretation 
as a cross-over phenomenon over a true, thermo- 
dynamically defined phase transition. 

3. CONCLUSION 

In the present analysis of the scaling of the 
maxima of the specific heat and susceptibility 
we found evidence for a cross-over regime in the 
four-dimensional Discrete Regge model of quan- 
tum gravity in the negative coupling region. If 
this can be substantiated by further investiga- 
tions and also for the Regge theory with contin- 
uous link lengths, the existence of a continuum 
limit at negative bare gravitational coupling is 
questionable. This is of major concern for a con- 
tinuum theory of quantum gravity with matter 
fields 
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